Model independent extraction of the pole and Breit-Wigner resonance parameters 



o 

(N 
4=> 

P-h 

in 



43 

I , 
Ph. 
<D ■ 
43 ' 



> 
OS 

O 

m 



x 

S3 



S. Ceci, 1 M. Korolija, 1 and B. Zauner 1 

1 Rudjer Boskovic Institute, Bijenicka 54, HR-10000 Zagreb, Croatia 
(Dated: February 18, 2013) 

We show that a slightly modified Breit-Wigner formula can successfully describe the total cross 
section even for the broad resonances, from light p(770) to the heavy Z boson. In addition to mass, 
width, and branching fraction, we include another resonance parameter that turns out to be directly 
related to the pole residue phase. The new formula has two mathematically equivalent forms: one 
with the pole, and the other with the Breit-Wigner parameters. 
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Resonances are unstable particles usually observed as 
bell-shaped structures in scattering cross sections of their 
decay products. For a simple narrow resonance, its 
fundamental properties correspond to the visible cross- 
section features: mass M is at the peak position, and 
decay width T is the width of the bell shape. These pa- 
rameters, along with the branching fraction x, are known 
as the Breit-Wigner parameters [l| . In reality, resonance 
peaks may be very broad, and the shape so deformed that 
it is not at all clear where exactly is the mass, or what 
would be the width of that resonance. In such cases reso- 
nance parameters are treated as energy dependent func- 
tions. These functions are often defined differently for 
different resonances. For example, in the case of p(770) 
resonance in the 7T7t channel, modern analyses include the 
pion-pion P-wave potential barrier (momentum to three 
halves) in the energy dependent width while in the 
case of Z boson the width function is proportional to the 
energy squared Q. 

With such model dependent parameterizations, the 
simple connection between physical properties of a res- 
onance and its model parameters is lost, and the choice 
of the "proper" resonance parameters becomes the mat- 
ter of preference. There are many definitions for Breit- 
Wigner mass, which is assumed by some to be the proper 
resonance physical property. Other will prefer the real 
part of pole position in the complex energy plane. Some 
will even define the resonance mass to be something un- 
related to these two most common definitions, as we will 
soon see, or assume that there is no difference between 
poles and Breit-Wigner parameters whatsoever. All that 
makes the comparison between cited resonance param- 
eters quite confusing and potentially hinders the direct 
comparison between microscopic theoretical predictions 
(such is and experimentally obtained resonance prop- 
erties ||. 

To clarify this situation we try to devise a simple 
model-independent formula for resonant scattering, with 
well defined resonance physical properties, which will be 
capable of successfully fitting the realistic data for broad 
resonances. 

In this letter we show how to dramatically improve the 



simple Breit-Wigner formula by incorporating in it just 
one additional (phase) parameter. This new formula has 
two equivalent forms that can be used to estimate cither 
pole or Breit-Wigner parameters in a model independent 
way. 

We begin our analysis by noting that resonant cross 
section is commonly parameterized by a simple Breit- 
Wigner formula [l[ 
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where q is a cm. momentum, J is the spin of the res- 
onance, while si and S2 are spins of the two incoming 
particles. Resonant scattering amplitude A is given by 



A = 



xT/2 



M - W-iV/2 



(2) 



where M is the resonant mass, T is the total decay width, 
x is the branching fraction to a particular channel (for 
inelastic scattering it is ^/x- m x out ), and W is cm. energy. 

This simple parameterization cannot describe most of 
the realistic cross sections since the resonance shapes are 
seldom symmetric. To fix this, a background contribu- 
tion is usually added. Unfortunately, there is no stan- 
dard way to add background, but polynomials in W 2 
(i.e., Mandelstam s) are commonly used in the literature 
(see e.g. @) 



\A\ 2 -> \A\ 2 + J2B k W 2k . 



(3) 
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To extract the resonance parameters, we do local fits 
of this parameterization to a broad range of data points 
in the vicinity of the resonance peak. To estimate the 
proper order n of polynomial background, we vary end- 
points of the data range and check the convergence of the 
physical fit parameters: M, T, and x. Goodness of the 
convergence is estimated by calculating c n j parameters 
for each data range and for all polynomial orders n and I 
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Smaller c n j means better convergence. 

To avoid false positive convergence signals as much as 
possible, we demand good convergence not just for two, 
but for three consecutive polynomial orders by using 



(5) 



Final result is the one having the smallest reduced \r 
among several fits (we use ten) with lowest convergence 
parameters c n . When statistical errors turn out to be 
unrealistically small due to the dataset issues, the spread 
in extracted pole parameter values is used to estimate 
parameters errors. 

To test this extraction approach, we analyze five broad 
resonances with well known properties, and masses rang- 
ing from less then 1 GeV, to almost 100 GeV. For 
A(1232) and AT(1440), we analyze GWU ttN elas- 
tic partial-wave amplitudes squared. For p(770) and 
Z boson we analyze e + e~ scattering ratio R (between 
hadronic an muonic channels) from PDG compilation Q , 
and for T(11020) the new BaBar data @]. 

Using the Breit-Wigner parameterization (|2J) on broad 
resonances does not produce very good results. There- 
fore, in the advanced approaches, resonance width V (and 
other parameters) are considered to be energy dependent, 
which drastically improves the fit. However, parameter- 
ization then becomes model dependent, obfuscating the 
connection between the model parameters and physical 
properties of the resonance. We want to find a simple 
model independent form, as close to the original Breit- 
Wigner parameterization as possible, that will be capable 
of successfully fitting the realistic data for broad reso- 
nances. To do so, we assume that the numerator and the 
denominator in Rel. ([2]) are functions of energy, expand 
them and keep only the linear terms. The new resonant 
amplitude A becomes 



A = 



M p -W-iT p /2 



\A B \e l 



(6) 



which turns out to be the lowest order Laurent ex- 
pansion of amplitude A about its pole position at 
W = M p — iTp/2. Therefore, M p and T p are pole mass 
and width, while x p T p /2 and 9 P are the complex residue 
magnitude and phase, respectively. (Note that we use the 
standard convention for residue phase 9 P , as used in PDG 
[j| , which differs from the mathematical residue phase by 
±7r.) Three additional fit parameters are residue phase 
6 P , (coherently added) background magnitude \Ag\, and 
the background phase 9b- We can extract only the rel- 
ative phase S p = 9 p — 9 B , since absolute square of this 
amplitude will be compared to the data. In order to ease 
the numerical analysis, we rewrite the new parameteri- 
zation in a compact form 



\Af 



\A f 



: (A*-W0 2 + A 2 

{Mp-Wy + Tl/t' 



(7) 



where fi and A are simple fit parameters related to the 
pole parameters through 



x p sin S p = \A B \ 
x p cos Sp = — \A 



iy2HAi 

T p /2 ' 
■ M p - n 
^~Tpj2~- 



(8) 
(9) 



Using parameterization ([7|1, we should be able to ex- 
tract pole mass, width, branching fraction, magnitude of 
the background amplitude, and the relative phase from 
the data. We again use the same polynomial background 
from relation ([3]) and convergence criteria from relations 
(@| and ([5]). However, at the very beginning of this anal- 
ysis we stumbled upon a problem with our fits. When 
we fitted A(1232) resonance, parameter A was rather un- 
stable, ranging from zero to several thousands MeV. In 
addition, fits often did not converge, even for a carefully 
chosen initial values. 

We looked into it more closely and realised that since 
A(1232) is almost elastic resonance (decaying by more 
than 99 percent to ttN channel), A should be zero due 
to the elastic two-body unitarity condition A^A = ImA. 
When A was set to zero, everything worked almost per- 
fectly. It is important to note that x p should be 1 for 
elastic resonances, again due to the unitarity, but setting 
A to zero does not imply that x p is 1. 

Things became really interesting when we tried to ex- 
tract Z boson parameters from e + e~ scattering data. The 
unstable A behavior seen in the case of A(1232) was ob- 
served again, even though Z boson is definitely not an 
elastic resonance. Fit could not be stabilized, and even- 
tually we tried A = again (x p still can take care of in- 
elasticity). This choice smoothed the fitting procedure, 
and the extracted resonance parameters were in excellent 
agreement with the PDG (pole) estimates [5[- Assuming 
that A = for other processes, we rewrite resonant am- 
plitude 
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with unmeasurable overall phase rj equal to 29b — 9 p . 
Square of this amplitude is then 



\Af 



[{M p - W) sinSp + Tp/2 cos S p Y 
(M p - W f + T2/4 



(11) 



We know that for A(1232), the overall phase ij iz zero 
due to unitarity, which means that 9b — S p and 9 P = 25 p . 
We compare our results for 25 p to published results of 
9 p for other analyzed resonances to check whether the 
same relation is valid for them as well. Roper resonance 
N(1440) is irN resonance with ttN branching fraction x 
estimated to 65%, and the 2S P value of —81° is a surpris- 
ingly close to the newest residue phase estimate —85° 
from Q. For the Z boson, these two values are even 
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closer: 28 p is —2.2°, while 8 p is —2.35°. Extracted masses 
and widths are much closer to the pole parameters listed 
in the literature, than to the Breit-Wigner ones, as can be 
seen in Table fl] The best fits for all analyzed resonances 
are shown in Figures [T][5] 

TABLE I: Resonance pole parameters extracted by using pole 
formula (|f f \ . Our 2<5 P is compared to the residue phase 9 P 
from the literature. PDG pole estimates are from Ref. [jj. 
The p meson pole and the Z boson residue phase are estimated 
by analytic continuation of Gounaris-Sakurai Q] and Breit- 
Wigner parameterization [f|, respectively. 



Resonance 


M p 1 MeV 


T p / MeV 


x p 1 % 


2S P 1 ° 


p(770) 


762 ± 1 


138 ± 2 


0.71 


1 ± 1 


POLE 


763 


144 


N/A 


N/A 


A(1232) 


1211 ± 1 


102 ± 1 


103 ± 1 


-47 ± 1 


PDG/POLE 


1210 ± 1 


100 ± 2 


104 ± 2 


-47 ± 1 


#(1440) 


1362 ± 5 


191 ± 10 


61 ± 4 


-81 ± 10 


PDG/POLE 


1365 ± 15 


190 ± 30 


65 ± 10 


-85 ± \l 


T(lf020) 


11000 ± 2 


43 ± 6 


0.10 


-52 ± 8 


BaBar [8J 


10996 ± 2 


37 ± 3 


N/A 


N/A 


Z(91188) 


91167 ± 6 


2493 ± 5 


15.4 


-2.2 ± 0.2 


PDG/POLE 


91162 ± 2 


2494 ± 2 


N/A 


-2.35 



Xb and 8b have almost exactly equal values as their pole 
counterparts x p and 8 P in Table fl] Furthermore, the ex- 
tracted pole and Breit-Wigner parameters are interre- 
lated through Manley relations [9( 

M b = M p - T p /2 tan 6 P , (15) 

r b = r p /cos 2 <v (16) 



TABLE II: Resonance parameters extracted by using new 
Breit-Wigner formula [|14[l. PDG estimates are from Ref. [jj. 



Resonance 


M b 1 MeV 


r 6 / MeV 


x b 1 % 


25 b 1 ° 


p(770) 


761 ± 1 


139 ± 2 


0.71 


± 1 


PDG 


775.5 ± 0.3 


146.2 ± 0.7 


0.69 


N/A 


A(1232) 


1233 ± 1 


120 ± 1 


102 ± 1 


-46 ± 1 


PDG/BW 


1232 ± 2 


117 ± 3 


100 


N/A 


AT(1440) 


1443 ± 2 


325 ± 11 


61 ± 4 


-80 ± 2 


PDG/BW 


1440 ± |8 


300 ± iBS 


65 ± 10 


N/A 


T(11020) 


11010 ± 2 


53 ± 8 


0.10 


-52 ± 8 


PDG 


11019 ± 8 


79 ± 16 


N/A 


N/A 


£(91188) 


91191 ± 5 


2494 ± 5 


15.4 


-2.2 ± 0.2 


PDG/BW 


91188 ± 2 


2495 ± 2 


15.3 


N/A 



The new pole parameterization (|11|) may be used in- 
stead of ([2]), since it works much better and adds only 
one quite important parameter S p . This parameter is the 
main ingredient of the shape of the resonance contribu- 
tion to the cross section. When 8 p is equal to zero, the 
new pole and the old simple Breit-Wigner parameteriza- 
tion are exactly the same. 

We began this study in the first place to find an 
improved Breit-Wigner parameterization, but ended up 
with pole parameterization instead. Following the origi- 
nal notion of Breit and Wigner, that the resonance mass 
is at the peak position and by noting the convenient 
form of our pole parameterization (jlpp that looks very 
similar to a single channel elastic amplitude (apart from 
x p 7^ 1 and 77 7^ 0). We now define the new Breit-Wigner 
parameters as single-channel K-matrix pole Mb, residue 
Fb, branching fraction Xb, and background phase 8b 



K = 



T b /2 
M h - W 



tan Sb, 



A = x b 



K 



1-iK' 



\A\< 



(r 6 /2 + tan5 6 



(M b - W) 2 + {T b /2 + tan 5b) 2 ' 



(12) 
(13) 
(14) 



In this form, Xb and 8b will be mathematically equal 
to x p and 8 P , respectively. When we fit parameteriza- 
tion (fl~4"l) to the data, extracted fit parameters Mb, Tb, 
and Xb are consistent with Breit-Wigner parameters in 
PDG [f| , as is clearly visible from Table flT] As expected, 
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FIG. 1: Fitting pole parameterization (|11[) to the data Q. 
The pole and BW masses have almost the same value, quite 
different from PDG estimate (dotted line). (We removed the 
data from the peak to eliminate the influence of the cj(782) 
resonance.) 

These Breit-Wigner parameters are uniquely defined 
and model independent, with directly observable mass as 
the peak of the squared amplitude \A\ 2 . However, they 
strongly depend on phase S p , which may change from re- 
action to reaction. That means that for the same pole 
position, there will be different Breit-Wigner masses and 
widths in different channels. Therefore, it is more prac- 
tical to have one (pole) mass and width for each reso- 
nance in particle data tables [j| than stockpiling differ- 
ent (Breit-Wigner) masses and widths for each process 
in which the resonance contributes. 
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FIG. 2: Fitting pole parameterization to the SAID data 
0| . Pole (solid) and BW mass (dashed) are clearly distinct, 
while the PDG estimates (dotted lines) are indistinguishable 
from them. All the data in this figure are analyzed, but only 
the black data points are used in the best fit. 
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FIG. 3: Fitting pole parameterization to the SAID data 
0|. This resonance has the largest difference between pole 
and BW mass. PDG estimates (dotted liness) are consistent 
with pole (solid) and BW (dashed) parameters. 
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FIG. 4: Fitting pole parameterization (fTTjl to the BaBar 
data Pole and BW masses are clearly distinct, but PDG 
estimate is neither of them. 
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FIG. 5: Fitting pole parameterization (fTTT) to the data [jj. 
The pole and BW parameters are consistent with their PDG 
estimates (dotted lines). 



In our study, there are two resonances that show sys- 
tematic discrepancy between the PDG mass estimates 
5] and our results presented here: the p(770) and the 
T(11020). For p meson an alternative parameteriza- 
tion by Gounaris and Sakurai [2] is used, where mass 
and width are defined somewhat unconventionally to 
take into account its mixing with w(782) and /9(1450). 
T(11020) mass and width were fit parameters to a Gaus- 
sian with a relativistic tail [lfj. For both resonances 
cited values in the PDG tables are neither masses con- 
sistent with the original Breit-Wigner idea, being at the 
peak of the resonance, nor the pole positions. Since the 
resonance parameters are collected in PDG tables to be 
used as an input for various models and for comparison 
between theory and experiment, placing all these reso- 
nance parameters in a single table may generally cre- 
ate considerable confusion. This confusion is evident in 
the p meson case where in the table with predominantly 
Gounaris- Sakurai masses (about 775 MeV) one can find 
pole masses (roughly 760 MeV). 

In conclusion, we have shown here that the original 
Breit-Wigner formula may be drastically improved by in- 
cluding a single additional (phase) parameter S p . Our re- 
sults suggest that parameter S p seems to be equal to the 
half of the resonance residue phase 6 P , regardless of the 
resonance inelasticity. This new formula has two equiv- 
alent forms that can be used to estimate either the pole, 
or the Breit-Wigner parameters in a model independent 
way. Having both forms enabled us to learn that in the 
PDG tables [5J there are values that do not correspond 
neither to pole, nor to Breit-Wigner parameters. Such 
an outcome undermines the proper matching between mi- 
croscopic theories (e.g., lattice QCD [4j) and experiment. 
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